function tech = technology(t,alpha, hc, I, PS, H, TFP_param)

TFP = exp(  TFP_param(1) + TFP_param(2)*(t) + alpha(4).*PS );

gamma1_tilde = alpha(6);
gamma2_tilde = alpha(7);

theta_tp1 =   ( TFP.*hc.^(alpha(1)).*H.^(alpha(2)).*I.^(alpha(3))  ).*(PS==1)  +  ( TFP.*[ gamma1_tilde.*hc.^(alpha(8)) +  (1-gamma1_tilde).*[ gamma2_tilde.*H.^(alpha(5)) + (1-gamma2_tilde).*I.^(alpha(5)) ].^((alpha(8))/alpha(5)) ].^(alpha(9)/alpha(8)) ).*(PS==0) ;


tech = theta_tp1;